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ABSTRACT 


An experimental investigation has been performed to establish 
the state of the gas in the nozzle expansion of air in a hypersonic 
arc-heated wind tunnel. Tests covered ranges of stagnation chamber 
pressure from k.20 to 9*3® atmospheres and stagnation enthalpies from 
700 to 3, ‘>10 Btu/lbm. Parameters measured during the tests were N 2 
vibrational temperature, N 2 rotational temperature, static pressure, 
static densi'v, velocity, total enthalpy, and hemisphere stagnation- 
point heating rate. Results of the te sts indicate general agreement 
with a chemical nonequilibrium, vibrational sudden-freeze model for 
the expansion of the air in the nozzle. 
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CHAPTER I 


INTRODUCTION 

The age of high-speed ballistic missiles and space travel has 
caused considerable interest in the phenomena occurring during hyper- 
velocity flight in planetary atmospheres, especially the earth's atmos- 
phere. The design of these vehicles requires a knowledge of the 
surrounding flow field. 

Although flight tests have been made to investigate the problems, 
ground facilities are used most extensively to acquire the needed data. 
Therefore, much effort has gone into developing wind tunnels that will 
simulate portions of the flight paths of various hypervelocity vehicles. 
Shock tubes, hotshot tunnels, expansion tubes, and arc-heated wind 
tunnels are among the devices that are used to study high-speed 
phenomena. 

A major problem encountered in simulating these high-speed flight 
conditions in wind tunnels is that thermal energy must be converted to 
kinetic energy by expanding the heated gas in a nozzle. Prior to 
entering the nozzle, the near-stagnant gas will be at a high temperature 
with various energy modes excited. For example, in the case of air, the 
total energy of the diatomic molecules will be distributed between the 
translational, rotational, and vibrational degrees of freedom, with 
energy also being absorbed in electronic excitation, dissociation, and 
ionization. As the gas expands down the nozzle, the expansion is 
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generally so rapid that the flow kinetics cannot keep pace with the 
rapidly changing thermodynamic state of the flow. Consequently, the 
gas dynamics of the flow is affected. 

Generally, translation and rotation will be represented by the 
same kinetic temperature due to the small number of collisions required 
to equilibrate these modes. However, vibrational and chemical changes 
require many more collisions to achieve equilibrium; therefore, chemical 
composition and vibrational temperature will generally not correspond 
to the equilibrium values that would be associated with the static 
temperature and pressure. In turn, the values of the flow variables 
(pressure, temperature, density, velocity, etc. ) are affected by the 
extent of nonequilibrium. 

These nonequilibrium phenomena present the researcher in the 
ground facility with two problems. First, the composition of the air 
in the free stream of the tunnel is generally quite different from the 
atmospheric free stream of a flight path. Secondly, the energy that is 
frozen in vibration and chemistry is not available for conversion to 
kinetic energy. 

Every effort must be made to determine the type of flow in the 
high-speed wind tunnels, and then the problem of nonequilibrium flow 
through shocks and over vehicles must be solved. These problems are 
rapidly becoming more severe as higher reentry velocities from distant 
planets must be considered. 

The present work is concerned with defining the type of flow in 
an arc-heated wind tunnel; that is, various experimentally determined 
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flow properties will be compared, with theoretical predictions from 
various gas flow models. Several authors have concerned themselves 
with this problem (Refs. 1 through 7). Theoretical methods for non- 
equilibrium flow in nozzles range from simple, sudden-freeze models to 
highly sophisticated models with coupling between chemistry and vibra- 
tion and chemistry and electronic excitation (Refs. 8 through 14). 

One can establish a priority of measurements that would contrib- 
ute to an understanding of the flow. They are, in decreasing order of 
relative importance and sensitivity, chemical composition, static 
temperature and vibrational temperature, static pressure, static 
density, and velocity. Presently, the determination of chemical com- 
position is ruled out since reliable spectroscopic techniques are not 
yet available for a flowing system (Ref. 7)* Ng rotational and vibra- 
tional temperatures can be measured using the electron beam Tow 
densities). Static pressures can be measured on the tunnel wall or on 
suitably designed probes. However, static pressure measurements can 
be highly uncertain in low-density flows. Density and velocity can be 
measured but are less sensitive to nonequilibrium effects. 

The measurements taken in the present study are N2 rotational 
and vibrational temperature, wall static pressures, static density, 
velocity, total enthalpy, and stagnation-point heating rate. The 
heating rates were taken in an attempt to lend credence to the values 
of total enthalpy used in the tests. 


CHAPTER II 


APPARATUS AND TESTS 
Tunnel 

The tunnel under investigation is the 1-Foot Hypersonic Arc 
Tunnel (Ref. 15) located at NASA’s Langley Research Center. Flow 
passes from the direct-current arc heater through a circular-arc throat 
into a 5° half-angle conical nozzle. At the 1- foot-diameter station of 
the nozzle, a cylindrical test section equipped with two quartz windows 
(1+ inches by 9*l/ 2 inches) is available for testing. The test section 
is followed by a 1-foot -diameter straight-pipe diffuser. A five-stage 
steam ejector provides the low back pressure for the expansion and 
exhausts the test gas to the atmosphere. The tunnel is a continuous- 
flow facility with test time limited primarily by electrode life. The 
tunnel is shown schematically in Figure 1. 

Power for the tunnel is supplied by a 900-kilowatt direct-current 
generator. An additional 900 kilowatts of power is available but was 
not used during the present tests. 

Some details of the tunnel components from the heater to the 
throat are presented to familiarize the reader with the specific com- 
ponents used, since there has been a progression of modifications to 
components after the publication of Reference 15. 

The cathc^e of the heater is a 1-3/^-inch-diameter doughnut- 
shaped inner electrode. This electrode is cast from a 20 per cent Cu 
to 80 per cent Ag alloy. Test gas is injected along the sting of this 
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electrode and through the arc. An orifice in the center of this elec- 
trode may be used either to measure pressure or to inject test gas. 

The outer electrode (anode) used with the tunnel depends on the 
enthalpy level desired. Two outer electrodes, D and E, are normally 
used. The E electrode has an internal diameter of 2.75 inches, while 
the D electrode has an internal diameter of 2.38 inches. Thus, E 
provides a larger arc gap with increased arc voltage and higher enthalpy. 
Both D and E electrodes have copper internal liners. The length of 
E electrode is 3 **l/ 2 in cbes and the length of D electrode is 3 inches. 

The plenum chamber that is used with E has a stainless steel 
liner with an inner diameter of 3*19 inches. The plenum chamber on 
D electrode is an integral part of the. electrode with an internal 
diameter of 3*25 inches. The plenum chamber length of D is 3 inches 
and that of E is 2-1/2 inches. 

Contamination traps are used with both outer electrodes to catch 
some (approximately 3 ° per cent) of the material eroded from the inner 
and outer electrodes (Ref. 16). The trap that is used with D has a 
1 -inch aperture, while the one used with E has a 1-3/^-inch aperture. 
The stagnation chamber pressure is measured with an orifice located on 
the downstream side of the trap. 

A gas injection ring is located just downstream of the contam- 
ination trap. Radially drilled holes equally spaced circumferentially 
inject the gas straight toward the center 01 the flow. This gas injec- 
cion ring is used to mix cooler gas with the hot gas that has passed 
through the arc. This procedure allows one to obtain a further 


7 


variation of enthalpy in addition to the levels produced hy the two 
different electrodes. The inner diameter of this ring is 2-1/2 inches 
and the length is 1-1/8 inches. 

The gas injection ring is followed by the circular-arc throat 
section. Throat diameters of 0.0624, 0.125, 0.142, and 0.2115 inch are 
available. The 0.2115-inch -diameter throat was used for the present 
tests. 

Two magnetic coils provide a 12,000-gauss axial magnetic field to 
spin the radial arc. The swirling of the arc insures electrode survival 
and stabilizes the arc. 

All tunnel components dowr to and including the throat section 
are water cooled. Inlet and outlet cooling water temperature and the 
water flow rate are measured on each component. 

Test Conditions 

Test conditions for the present experiments will be given in 
terms of input parameters to the tunnel and the resulting output tunnel 
test conditions. 

The total air mass flow through the tunnel was held constant at 
0.024 lbm/sec. From 0 to 50 per cent of this flow was room-temperature 
dry air added through the injection ring when a variation in enthalpy 
was desired. The arc current was held at 1,500 amperes when outer elec- 
trode D was used. Electrode E was operated with 1,200 amperes 


arc current. 
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With these inputs, stagnation chamber pressure varied from 
4.20 to 9.38 atmospheres and stagnation enthalpy varied from 700 to 
3,510 Btu/lbm, thus giving a stagnation chamber temperature range from 
2,665° to 8,225° R* The test section parameters of Mach number, 
Reynolds number, etc., are the objects of the tests and will be dis- 
cussed later in the text. However, pitot pressure at a location 
2-1/4 inches downstream of the nozzle-test section juncture varied 
from 6.3 X 10“3 to 1.2 X 10“^ atmospheres. 

Models and Instrumentation 

The models and measuring devices used during the present tests 
are listed below. 

(1) A total calorimeter for the determination of tunnel total 
enthalpy. 

(2) The electron -beam apparatus for the measurement of N 2 
rotational and vibrational temperatures. 

(3) Static pressure orifices in the nozzle wall. 

(4) A mass flow probe for the measurement of free-stream 
density and velocity. 

(5) A 1-inch nose -radius — hemispherical stagnation point heat- 
transfer probe with a 1-inch-diameter cylindrical afterbody. 

(6) Varxous pitot probes in conjunction with (2), (3), (4), 
and (5) above. 

Total calorimeter .- A total calorimeter was used to measure 
total enthalpy because delays in the delivery of a free-stream total 
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enthalpy probe prevented its use in time for inclusion in this paper 
and because an energy balance on the arc heater was too inaccurate a 
method for this facility, since the heater efficiency was less than 
25 per cent. The calorimeter was attached to the downstream end of 
the tunnel throat section and the first two nozzle sections of the 
tunnel downstream of the throat were removed during the calorimeter 
tests. 

The calorimeter was constructed of coiled copper tubing sealed 
with soft solder. The calorimeter inner diameter decreased from an 
8-inch-long, 3“ihch-inside -diameter entrance section to a constant 
1-1/2 -inch-inside-diameter section 12 inches long. An insert apparatus 
was used to decrease the flow velocity and to make more surface area 
available for heat transfer from the gas. Gas residual temperature was 
measured at the exit of the calorimeter with stainless- steel- sheathed 
chromel-alumel and iron-constantan thermocouples placed at various 
radial positions. The signal from these thermocouples was recorded on 
a 52-channel dry-develop oscillograph. 

The inlet and outlet temperatures of the calorimeter cooling 
water were measured with iron-constantan thermocouples, and these 
temperatures were read visually by balancing a potentiometer. Calorim- 
eter cooling water flow rate was measured with a turbine -type flow 
meter and the output was recorded on a recording potentiometer. 

Electron beam apparatus .- A 15-kilovolt electron-beam gun and a 
l/2-meter scanning spectrometer were used to measure rotational and 
vibrational temperature of N2 in the test section 13 inches from the 
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nozzle exit. A 1-inch-diameter hemisphere-cylinder pitot probe was 
located l/k inch downstream from the beam during this series of tests. 

The chamber of the gun was maintained at a pressure of approxi- 
mately 10' 6 torr by a liquid nitrogen-trapped diffusion pump. The beam 
width was approximately 0.060 inch and was passed in the top of the 
tunnel and across the stream (Fig. 2). The tunnel suited as the electron 
catcher. The beam current was maintained constant at 200 microaraps. 
References 17 and l8 describe the electron -beam gun used in these tests. 

A l/2-meter scanning spectrometer was used to view the beam. 
Entrance and exit slit widths were adjusted to resolve the rotational 
lines of the (0,0) vibrational band in making the rotational temper- 
ature measurements. These slits were opened to get a gain in intensity 
so that the (0,1) Ng and (1,2) Ng vibrational bands were resolved in 

making the vibrational temperature measurements. A grating with 

o 

2,400 lines/mm (resolution of 0.2 A) was used on the rotational temper- 
ature studies and a grating with 600 lines/mm (resolution of 0.8 A) was 
used for the vibrational temperature studies. 

Wall static pressure measurements . - The two static pressure 
orifices were located 6 and 6-1/2 inches upstream of the nozzle-test 
section juncture, respectively. One orifice was located directly above 
the tip of the pitot probe and the other was located 5° upstream of 
this point to check source flow effects. Location of the orifice up- 
stream of the nozzle-test section juncture was required to minimize 
the effects of any disturbance from this juncture. The pitot probe 
used for the static pressure tests was a l/4-inch-diameter flat-faced 


0 to 50 kV 
power supply 


-Plexiglass 


Control grid - 
Cathode - 

Focusing anode- 
Accelerating anode - 


* Vacuum 


-Iron $hi*k$ 


Ammeter 

and 

recorder 


Spectrograph _/ / * 

Quartz window^ 

Electron beam. 

Tunnel test section. 


Figure 2.- Schematic drawing of the electron -beam spectrometer 
apparatus attached to the tunnel. 
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cylindrical probe 12 inches long. The probe was sufficiently small to 
avoid any effects of probe size on the static pressure measurements. 

Static pressures and all pitot pressures were measured with 
capacitance -type transducers. Sensing heads of 10, 5 0 > and 100 torr 
were used. The signal conditioners for these gauges provide ranges 
of 1.0, 0.3> 0.1, 0.03, 0.01, 0.003, 0.001, and 0.0001 of the stated 
head pressure. These gauges are very accurate, calibrations indicating 
better than 2 per cent (of the readings) agreement with the calibration 
standard (McLeod gauge). Comparison of the gauges on a common manifold 
indicated agreement between the various gauges ranging from 5 per cent 
at 0.030 torr to much better than 1 per cent at higher pressures. 

When not in use, these gauges were at all times maintained at 
approximately 10"^ torr. Output from the pressure gauges was recorded 
on the 52-channel oscillograph. 

Mass flow probe .- The shock- swallowing mass flow probe is shown 
in Figure 3. The design of this water-cooled probe is based on the 
results of the studies reported in Reference 19* Equal interior and 
exterior tip angles of 30° were used to insure that the stream tube 
which was swallowed was equal in cross-sectional area to the geometric 
inlet area of the probe (inlet diameter = 0.240 inch). The inlet 
section was followed by a l/4 -inch-long constant area section with a 
length-to -diameter ratio of 1.22. The constant area section was followed 
by an expanding region to a diameter of 0.431 inch. The remainder of 
the probe was of this diameter. This design of the entrance section 
gave the best results of several tips tested in Reference 19. The tip 
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design is of great importance, since viscous or bluntness effects can 
change the effective entrance area. 

Two copper-constantan thermocouples were located near* the exit 
of the 0.1+30 inner diameter tube. One thermocouple measured wall 
temperature, while the other measured the centerline gas temperature 
of the flow through the probe. The air passed the thermocouples, 
flowed through a piece of 1-inch I.D. flexible tygon tubing, through 
a flow rate measuring orifice, and back to the low static pressure of 
the tunnel. The tygon tubing allowed the probe to be moved so that 
lateral surveys could be made. 

A 0.187-inch -diameter corner-tap orifice was used to measure the 
probe mass flow rate. The orifice was calibrated in the probe with the 
probe installed in the tunnel test section. Dry air was used for the 
calibration and a calibrated rotameter was used as the flow rate 
standard. Supersonic flow across the orifice was verified for the 
expected range of probe mass flow rates. Also, measurements of pres- 
sure drop through the probe during the orifice calibration indicated 
that the shock would probably be swallowed during the tests. 

The capacitance pressure gauges discussed previously were used 
to measure pressure upstream and downstream of the orifice. These 
gauges were also used to measure pitot pressure when a valve downstream 
of the orifice meter was closed. 

The 1-inch outer diameter of this water-cooled probe is com- 
patible with the probe size used during the rotational and vibrational 
temperature studies. The mass flow probe was located at the same 
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position in the tunnel as the electron beam (13 inches from the nozzle - 
test section juncture). 

The inlet orifice diameter of the mass flow probe was equal to 
the tip diameter. Thus, bleed-off effects (lowering of the pressure) 
when the probe was used as a pitot probe were a possible error source. 
Another pitot probe with the same nose shape, but with an orifice - 
diameter-to-t ip-diameter ratio of 0.333# was built so that these pitot 
pressure measurements could be compared with those from the mass flow 
probe . 

Stagnation point heat-transfer probe .- The stagnation point heat- 
transfer measurements were taken on a blunt cylindrical probe with 
R n /R b = 2.0 and R c /R b = 1.0. The probe tip was located 2-1/4 inches 
from the exit of the nozzle -test section juncture. 

The principle of operation of the probe was that of the Gardon 
gauge (Ref. 20). A 0.130-inch-diameter const antan disk with a thick- 
ness of 0.015 inch was soldered to the surface of a cavity in the 
tip of the blunt probe (see Fig. 4). A 30-gauge copper wire was 
attached to the center of the disk. The cavity behind the constantan 
disk was a sealed insulating air pocket. Thus, cooling of the disk 
had to be accomplished by radial heat flow to the copper heat sink. A 
differential temperature exists from the center of the disk to the 
circle of contact with the copper when heat is applied to the tip of 
the probe. Ihus, the heat-transfer rate is proportional to this AT. 
This heat-transfer probe was calibrated on a radiation calibration unit. 
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A 1- inch-diameter pitot probe exactly like the heat-transfer 
probe was used on the heat-transfer tests. Hie probe was located at 
the same position in the test section as the heat-transfer probe by 
alternating insertion of the two probes into the stream. 

Test Procedure 

Prior to a test, the particular probe in use was carefully 
aligned with the flow on the tunnel centerline. In the case of the 
mass flow probe, the tip region was thoroughly cleaned before each test. 
Calibrate signals were applied to all associated instrumentation to 
insure that the equipment was functioning and that the calibrations 
had not changed. 

The tunnel was evacuated to approximately 0.020 torr prior to 
a test. Cold flow was then passed through the tunnel to determine if 
there were any air leaks in the heater section. Hie arc-starter tab 
was also checked at this time to verify electrical contact between the 
inner and outer electrodes. With the cold mass flow passing through 
the tunnel, the arc was fired and the arc current brought up to the 
desired level. Tunnel-associated measurements (P^. etc. ) required no 
more than 10 seconds to stabilize. Test times ranged from approximately 
30 seconds on the heat transfer and calorimeter tests to approximately 
4 minutes for a rotational temperature measurement. 

The calorimeter tests consisted of allowing the air-out and 
water-out temperatures to come to equilibrium. Hien, the potentiometer 
measuring the water-out temperature was balanced and read visually. 
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The vibrational and rotational temperature measurements involved 
scanning the desired range of wave lengths. More than one scan was 
taken if any data irregularities were noticed. Pitot pressure was 
measured on the centerline near the end of the test. 

Wall static pressures were generally measured with the i/4-inch 
diameter pitot probe located on the tunnel centerline. Visual monitoring 
of the pressure (by use of the visual dial on the signal conultioner) v. t s 
required to determine when the static pressures had stabilized. 

The mass flow probe was operated with the valve open and flow 
passing through the probe until the upstream and downstream pressures 
on the flow-rate measuring orifice had stabilized. Then, the valve was 
closed and pitot pressure was measured. 

For the stagnation point heat-transfer tests, the pitot probe was 
first inserted into the stream and pitot pressure was measured. Ibis 
probe was th~n withdrawn and a measurement of heat transfer was taken 
by inserting the stagnation point heat-transfer probe. 

After a test was completed, all capacitance pressure gauges were 
again switched to the vacuum pump. A cold mass flow equal to the hot 
mass flow was passed through the tunnel and cold stagnation chamber 
pressure was recorded for use in the sonic-throat enthalpy technique. 


CHAPTER III 


DATA REDUCTION 


Calorimeter 

The purpose of performing tests with the total calorimeter was to 
determine the total enthalpy of the flow at the entrance to the expan- 
sion (i.e., at the entrance to the convergent throat section). Any 
difference between the total enthalpy at this point and in the test 
section of the tunnel depends on the validity of the assumption of an 
adiabatic core of flow down the tunnel. An energy balance on the throat 
section, the calorimeter, and the air out of the calorimeter gives the 
total enthalpy at the entrance to the expansion as 


H t,CAL = H t,l 


(i "°Pw ^throat + (ilwC Pw At) CAL 




CpgTa, out 


where the subscript a refers to air and the subscript w refers to 
cooling water. 

The equilibrium sonic-throat technique for calculating total 
enthalpy was also used so that a correlation between this method and 
the calorimeter results could be obtained for use in later tests. 
Equations for total enthalpy by the sonic-throat method in equilibrium 
air are given in Reference 21 as functions of tunnel mass flow rate, 
throat area, and plenum chamber pressure. For the present tests, graphs 
were used to determine enthalpy at the measured pressure as a function 


19 




20 


p. . 

of .-x i . S -— /T + . (Fig. 5)* Thermodynamic data used for this plot 

p t,c m h 

were taken from the equilibrium throat conditions of the Cornell Aero- 
nautical Laboratory computer program (Ref. 12). Similar plots are 

presented in References 21 and 22, where H^. is plotted against 
P A* ’ 


A , in practical usage, must be the effective throat area and 

^ p + ^ m p. ,a* 

the use of — £ J T + _ in place of — — indicates that the 

r t, c m 

assumption has been made that A*^ c = and the validity of the 

sonic-throat technique is dependent on this assumption. 


Vibrational Temperature 

The electron beam technique was used to measure free- stream 
nitrogen vibrational temperature in air. The vibrational temperature 
was determined by comparing the relative intensities of the (0,1) and 
(1,2) vibrational bands of A plot of the relative intensity ratios, 
of these bands is shown in Figure 6 as a function of Ng vibra- 

U,2> 

tional temperature (Ref. 17). The fact that these two bands are 
reasonable to use has been confirmed by Petrie (Ref. 5) and others 
(Refs. 17 and 23). 

A typical trace obtained from the scanning spectrometer is shown 
in Figure 7. The ratio of the intensities of the vibrational bands was 
determined by ratioing the areas contained within the bands. If any 
irregularities occurred during the test (such as beam fluctuations), 
the test was prolonged and the scan was repeated. 
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Rotational Temperature 

The relative intensities of the rotational lines of the R-branch 
of the (0,0) vibrational band of N* were used in the determination of 
the rotational temperature. A typical spectrometer scan of the rota- 
tional lines within this band is shown in Figure 8. 

The most commonly used model for calculating rotational tempera- 
ture from the measured relative intensities is that of Muntz (Ref. 2k). 
The intensities of the rotational lines according to this excitation- 
de-excitation model can be expressed as 



K" 


K',K" 


- k '( k ' + 1 ) V/Te 


( 2 ) 


The rotational temperature can be obtained by plotting 

in [ I K« # K"/ K,Gvlf K',K" ^ a ® ainst K '( K ' +1 )‘ v k',k" iS the wave 
number and is very nearly constant for the band. G is a function of 

K’ and T^ and therefore the solution far T^ is iterative. The slope 

of the straight line on the intensity plot is where 0^^ 

for N 2 is 2.891° K. 

Several problems are often encountered when attempting to reduce 
the intensity data from the spectrometer. Ashkenas (Ref. 25) presents 
a summary of these problems along with a critical review of data reduc- 
tion techniques. Seme of these problems are: 

(1) Overlap of the P-branch onto the R-branch of the (0,0) band 
of N+ with possible incorrect intensity values. 

(2) Poor resolution of rotational lines resulting in incorrect 


intensity values. 
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(3) Overlap of other bands, limiting the number of useful 
rotational lines. 

(4) Beam current and/or density variations during the spectrum 
scan, resulting in erroneous intensity values. 

(5) Failure of the measured intensity values to define a 
straight line when plotted, as mentioned above. 

As noted by Ashkenas, (l) and (3) pose no serious threat to the 
method; (2) is a function of signal strength and type of spectrometer 
and is not a serious problem at higher rotational quantum numbers; 

(4) may be avoided by careful monitoring during the test; ( 5 ) non- 
linearity of the data poses a limitation on the technique and does not 
seem to have been resolved satisfactorily in the literature. 

In most research investigations, the higher rotational levels 
seem to be overpopulated, thus producing the nonlinear effect. The 
possibility of not being able to define a rotational temperature then 
arises. According to Reference 25, if the gas is in thermodynamic 
equilibrium with no overlapping, all measured spectral lines must fall 
on a straight line on the intensity plot. If this does not occur, then 
the excitation-de-excitation model must be questioned. 

Another problem that arises with the Muntz model is the consis- 
tent overprediction of the true temperature; that is, the predicted 
values are higher than are obtained from other methods (thermocouples, 
inference from pitot pressure in flowing systems, etc. ). This problem 
becomes much worse at lower temperatures. 
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As mentioned previously, both of these problems are discussed 
in the literature, and various solutions are proposed and reasons 
hypothesized. Robben and Talbot (Ref. 26) encountered both problems 
(i.e., nonlinearity of the intensity curve and the overprediction of 
rotational temperature by the Muntz model). These authors suspected 
that secondary electrons could be the cause of the problem. They pre- 
sented a rotational temperature correction curve based on their data 
obtained in nozzles and free jets. 

Marrone (Ref. 27), based on his own free -jet experiments, con- 
tended that the high temperatures encountered by Robben and Talbot 
(Ref. 26) could be due to rotational freezing, although he, too, 
recognized secondary electrons as a possible cause. He also encountered 
the nonlinearity problem and questioned whether it might be due to 
rotational freezing or to secondary electrons. Tirumalesa (Ref. 28) 
concluded that the high temperatures can be explained by rotational 
freezing. 

Hickman (Ref. 29) suggested that the model proposed by Muntz 
was incorrect, and he presented a new model for the reduction of rota- 
tional temperature data. This new model allows rotational changes of 
AK = tl and +3 as opposed to the Muntz model which only allows 
AK = tl. According to Hickman, his model will give accuracy to within 
the magnitude of the characteristic rotational temperature of N 2 « Data 
reduction according to Hickman's model also utilizes Equation (2), but 
with a different expression for G. 
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Ashkenas (Ref. 25) compared Muntz's and Hickman's models with 
thermocouple data obtained in a near- static gas at 144° R and 540° R 
for a range of densities. He reports similar shortcomings of both 
models, although Hickman's model does give lower temperatures. 

Ashkenas found that a rotational temperature measurement is dependent 
on the density and on the number of spectral peaks used in an analysis. 
Limited corrections (near 144° and 540° R) could be made to rotational 
temperature data on the basis of Ashkenas' work; however, some degree 
of scatter is evident in the correlations that would have to be used. 

Lillicrap (Ref. 50) noted that the model used by Muntz predicts 
rotational temperatures that are too high and suggested secondary 
electrons as the cause since they will not necessarily obey the optical 
selection rules. He successfully reduced the data of Ashkenas (Ref. 25) 
with a formulation including the effects of secondary electrons. How- 
ever, the application as presented is essentially limited to the data 
of Ashkenas. 

Because of the uncertainties involved in obtaining rotational 
temper ature from the intensity data, the experimental data was reduced 
using both models. However, on the basis of Ashkenas' work (Ref. 25), 
Hickman's model should give the best results in the density range of 
the present tests. 

In reducing the data, any spectral line intensities that were 
suspected of being affected by poor resolution, overlap of the P-branch, 
or electronic noise were omitted. This generally meant not using the 
first two and the last two spectral lines appearing in the band. The 
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scatter and nonlinearity of the remaining points was reasonable and a 
straight-line fit was made to this data on an intensity plot. 

A typical rotational temperature data reduction plot is shown 
in Figure 9* These data points are from the spectrometer trace shown 
in Figure 8. The faired line is from a computer program modeled after 
that of Hickman (Ref. 29) which was used to check the data reduction 
results. This program uses a least squares fit to the data and iterates 
seven times to insure convergence of the iteration process. 

Wall Static Pressure 

A nozzle distribution of pressure measurements was not taken due 

, c 

to the 5“l/2-foot length of the shallow 5 half -angle nozzle. Two 
static pressure orifices were used; one was located directly above the 
pitot probe tip and the other was located ahead of the tip as a check 
on source flow effects. Readings from these two orifices generally 
agreed within 3 per cent or less. 

Many investigators have used slender ogive-cylinder probes with 
orifices located well back from the tip to determine free-stream static 
pressure (Refs. 1, 2, 3, and k). This technique involves a calculation 
of the effects of boundary-layer-induced pressure (Ref. 31) in the 
measured pressure. 

Also, in high-temperature, low-density flows, thermal transpir- 
ation (Ref. 32) and the heat-transfer orifice effect (Refs. 33 and 3*0 
may influence the pressure measurement. All of these factors must be 
considered for tests in the 1-foot hypersonic arc tunnel. 
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Several difficulties also arise when using nozzle wall pressures 
as a diagnostic technique. One problem is concerned with taking meas- 
urements in the low-density regions. As the boundary layer along the 
nozzle grows, a contouring effect occurs. This is illustrated in 
Figure 10, where the equilibrium boundary- layer growth for three 
different enthalpies is shown. These theoretical data were taken from 
an existing program. The wall static pressures may be increased in the 
low-density regions by weak compression waves created by the boundary- 
layer contouring. Petrie (Ref. 5) and others have noted a relative 
increase in wall static pressure as the low-density regions were reached. 
MacDermott and Marshall (Ref. 7) measured wall static pressures that 
were higher than theory and corrected their data by equating wall 
pressure to pressure at the edge of the boundary layer and relating 
the pressure at this location to that at the nozzle centerline by use 
of measured pitot profiles and assumed isentropic flow in the core. 

No such corrections were made to the present data, since such a correc- 
tion technique requires a better definition of core size than was 
apparent from the pitot surveys (Fig. 11). 

A correction was made to the pressure data for the heat- 
transfer orifice effect. This phenomenon involves an erroneous low 
pressure indication due to orifice size in a low-density flow with 
heat transfer. In th '• present tests, upward correctiors to the meas- 
ured data ranged from 3 per cent at the lowest enthalpy to 12 per cent 
at the highest enthalpy. Thermal transpiration effects were not 
important since the wall temperature was the same as the gauge temperature. 





Distance from centerline, inches 


Figure 11.- Pitot pressure surveys across the test section for arrange 
of test conditions (surveys taken 2-l/b inches from nozzle exit). 



Pitot Pressure 




Pitot pressure data were readily reduced except for a second- 
order viscous effect in the merged layer regime. Corrections for the 
viscous effect on pitot pressure found in References 35 and 36 indi- 
cated that this correction was less than 2 per cent and therefore should 
be neglected. Agreement between probes ranging from l/4-inch diameter 
to 1-inch diameter also supported the neglecting of the viscous 
correction. 

Density and Velocity 

The data obtained from the mass flow probe were probe mass flow 
rate (with the valve open) and pitot pressure (with the valve closed). 
The probe mass flow rate is equal to the product of free- stream density 
and velocity and the entrance area of the probe, that is. 


nip PooVcoAp ( 5 ) 

if the shock is swallowed so that the cross-sectional area of the 
swallowed stream tube is equal to the geometric area of the probe 
entrance. The equal interior and exterior angles of the tip were 
designed to achieve this. It is difficult to prove that the probe 
swallowed the shock; however, photographs of the probe tip with the 
valve closed and with the valve open (Fig. 12) appear to indicate that 
the shock is swallowed. The self-luminescence of the flow made the 
shock visible. Note the presence of the gas cap with the valve closed 



(a) Valve open. 



(b) Valve closed. 

Figure 12.- Photographs of the mass flow probe tip illustrating the 
disappearance of the stagnation region gas cap when applying a 
very low downstream pressure to the probe. 







36 


and its absence when the valve is open. When the valve was slowly 
opened, the slow recession of the gas cap into the probe could be 
observed. The only other proof of the swallowing of the shock is an 
indirect proof by the comparison of measured density and velocity with 
theory. 

No corrections were necessary to the value of probe mass flow 
rate obtained from the sonic orifice. Probe centerline gas tempera- 
tures upstream of the orifice indicated that the flow was cooled to 
the same temperature at which the orifice was calibrated. Thus, with 
a knowledge of the probe flow rate and the entrance area (d = 0.2^0 in. ), 
the product of free-stream density and velocity was obtained from 
Equation (3)* 

In hypersonic flow, the pitot pressure may be closely approxi- 
mated by the product of free-stream density and free-stream velocity 
squared, that is, 


P t ,2 = BP.vi (4) 

The factor B is very weakly dependent on the ratio of specific heats 
and, consequently, is weakly dependent on the particular gas model 
chosen (equilibrium, nonequilibrium, or frozen). The equation 

P t,2 - °-9 £ 5 P jl 


(5) 
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was used since the approximate gas model was known from the rotational 
temperature measurements. 

Comparisons of the pitot measurements from the mass flow with 
the pitot measurements from the geometrically similar prohe with 
d/D = O.553 indicated as much as a 3 per cent upward collection to 
the pitot pressure obtained from the mass flow probe. Thus, the pitot 
pressure measured by the mass flow probe was relatively unaffected by 
pressure bleed-off effects. 

Equations (3) and (5) may be combined to obtain the expression 
for free- stream velocity as follows: 


V« = 


ts. 

0.925 


■JbL 


( 6 ) 


Free -stream density may be calculated by inserting this value of V,*, 
into Equation ( 3 ). 


Stagnation Point Heat Transfer 

In order to reduce the heat-transfer data, an effective radius 
of the hemisphere was required since the body radius was 0.5 inch and 
the nose radius was 1 inch (not a true hemisphere). The probe was 
originally intended to have a 1 -inch-diameter hemispherical tip; how- 
ever, installation of the constantan disk blunted the tip. The true 
shape of the probe was determined by a 10 power magnification on an 
optical comparator. The actual shape of the probe is shown in 
Figure 13 by the solid lines. 
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References 37, 38 , and 39 present data on st agnation -point 
velocity gradients for nose shapes varying from hemispheres to flat- 
faced probes. Using these velocity gradients, a particular shape can 

a 

1 

be simulated by a hemisphere with an effective radius. eff is a 
function of nose radius, body radius, and corner radius. For all of 

i 

these studies, the center of the corner radius was located on the radius 
arm of the nose radius at the intersection point of the two arcs; that 
is, the corner arc and the nose arc had the same slope at the point of 
intersection. 

A region of rapidly changing curvature formed the transition 5 

from the nose radius of curvature to the comer radius of curvature in | 

the present case. To this author's knowledge, no information on 
stagnation-point velocity gradients on this type of configuration is 
available. Thus, an effort was made to form reasonable estimates of 
the limits of the effective radius from the data presented in 

Reference 39 * : 

4 

The first obvious limit on %^ e ff is that it can be no less J 

than the radius of the original hemisphere, that is, %^ e ff > i/ 2 i nch ' '! 

The other l imi t is more difficult to determine. The effective corner f 

radius which would be tangent to the nose axe at the intersection with 
the original comer radius was determined graphically. The other con- 1 

' I 

straint on this effective comer radius was that it become tangent to | 

the 1-inch-diameter cylindrical afterbody extended. The effective | 

corner radius was found to be O.359 inch. 

I 
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This new body now was of the same class studied in Reference 39* 
Thus, another limit for the effective hemispherical radius was found 
from this reference as a function of nose radius, effective corner 
radius, and body diameter. %,eff determined to be 0*740 inch. 

Thus, 

0.500 < < 0.?40 

are the limits (in inches) on the effective hemispherical nose radius. 
These are still rather broad limits and restrict the ability to com- 
pare the heat-transfer data with theory. However, the stagnation point 
heat -transfer rate is inversely proportional to the square root of the 
hemisphere radius and, therefore, the uncertainty in heat-transfer rate 
is considerably less than the uncertainty in the radius. 


CHAPTER IV 


THEORY 

Enthalpy Determination 

Plenum chamber .- The total enthalpy at the exit of the plenum 
chamber was measured using the total calorimeter. The enthalpy at 
this point was also determined theoretically by use of the equilibrium 
sonic-throat method. The validity of the sonic-throat method depends 
upon a knowledge of the thermochemical state of the gas in flowing from 
the plenum chamber to the throat as well as a knowledge of the effec- 
tive throat area. The theory of the sonic-throat technique, including 
some factors which may affect the method, is covered extensively in 
Reference 21. Reference 22 also provides good coverage of the method. 

Only the equilibrium sonic-throat method will be considered here, 
since the technique was used as a correlation tool in the present tests. 
The equations for the steady, isentropic, one -dimensional flow of a gas 
through the nozzle are 

Continuity m = PAV (7) 

Energy = h + ^ (6) 

State P = ZpR g T (9) 

At the nozzle throat, 


m = p*A*V* 
4l 


( 10 ) 




k2 


Combining Equations (8), (9), and (10), one obtains 


m 




A*P. 


t,l 


z t,i R s T t,A p t,i 



21 - 


«t,l, 


(ID 


For a perfect gas, this reduces to (Z^. #1 * 1): 


_J> i_ Cl/^_V + 1/r ' 1 

A * p t,i v ’-H* + x > 


( 12 ) 


and the dependence of total enthalpy on mass flow rate, effective area 
ratio, and plenum chamber pressure is noted. 

Winovich (Ref. 21) derives expressions from the general equation, 
(11), for a perfect gas, a c&lorically perfect gas, and a real gas. 

Po. HI 

However, in the present case, plots of — ■ *- yT^ c versus 

t, c 

were used for the sonic-throat method (Fig. 5 )* The values of 
p. . m . 

■ ■X ? 1 c were obtained from the equilibrium throat conditions 

P t,c m li 

of the Cornell Aeronautical Laboratory computer program (Ref. 12) by 
using the assumption that A*f f ^ c = A*f^ h . 

Test section .- The enthalpy in the test section of the wind 
tunnel is the same as that in the plenum chamber if the gas has expanded 
adiabatically through the nozzle. This assumption is commonly made, 
although it should be checked in an arc-heated wind tunnel. Probably 
the best technique for determining the total entnalpy in the test 
section is the shock- swallowing enthalpy probe with separate cooling 
passages for the internal and external coolant flow. However, in the 
absence of such a probe, the stagnation point heating rate is often 
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used to infer the enthalpy. This is accomplished by relying on a 
theoretical or empirical relation involving the heat-transfer rate and 
the total enthalpy. For example, Zoby (Ref. 40 ) presents the following 
empirical formula for stagnation point heat transfer to a hemisphere in 
air. 


<1 



0.046l(H. . 

Zyl 



(13) 


where q is in Btu/ft 2 -sec, R H is in feet, p t,2 is ^ a^ospheres, 
and ^ and h^. are in Btu/lbm. This equation inherently assumes 
a fully catalytic wall with a Lewis number of unity. 

In low-density, high-enthalpy flows, one must exercise extreme 
care in using the stagnation point heat-transfer rate to determine 
total enthalpy. Two problems which affect the measured heat-transfer 
rates are the catalytic efficiency of the gauge surface and second- 
order viscous effects. If the stagnation point boundary layer were in 
equilibrium and viscous effects small, the data would still be useful 
for calculating total enthalpy. In any event, it is appropriate to 
determine what conclusions can be drawn from the data, even with the 
uncertainty in effective radius. 

Rosner (Ref. 4l) presents an equation for stagnation point heat 
transfer to a hemisphere, including the effects of varying degree of 
wall catalycity. This equation is identical to the frozen boundary- 
layer equation of Fay and Riddell (Ref. 42) with the exception of the 
term involving the wall catalycity. This equation is 


4 . 0.763^) (PP e ^)°- 5 (Pr)-°^ ♦ 0|(Le) 0 - 6 - 


where 


^ ^t,l " ^f,w^ ^ 0^chem,max 


(lb] 

(15 




(for s. completely frozen 
■boundary layer) 


S = 0.763(3P e ^ e )°* 5 (Sc)- 0 - 6 p/ 1 


(16 


(17 



In the preceding equations, the extent of recombination, 0, 
contributes perhaps the most uncertainty to the calculated value of 
heat transfer. If both the shock layer and the boundary layer are in 
a state of nonequilibrium, the calculation of 0 becomes a monumental 
task. If conditions at the edge of the boundary layer are known and 
if the boundary layer is frozen, 0 may be expressed in terms of the 
surface catalytic reaction rate constant, K^, and the quantity S, 
which is a measure of the diffusion velocity in the boundary layer. 
Even in this case, Ky is a relatively uncertain factor. for 


cupric oxide was calculated for use with the present tests, using data 
from Reference 43 ; however, extreme caution must be exercised in using 
values in the literature (Carden, Ref. 44) since surface catalytic 


V V “ 


properties of metals are dependent on surface preparation as well as 
the duration of exposure to atomic specie. 

Prior to applying the heat-transfer theory, some knowledge of 
the state of the gas in the shock layer and in the stagnation point 
boundary layer is necessary (i.e., Is the gas frozen? In a nonequilib- 
rium state? Or in equilibrium? ) This question may be investigated 
by calculating the Damkohler number which is defined as 


Da 


Residence time 
Reaction time 


( 19 ) 


In the case of the shock layer, the residence time of the fluid 
is taken to be 2 Rjj/Voo. Reference 46 presents the following expression 
for the Damkohler number for chemical dissociation in the shock layer: 



( 20 ) 

where T is °K, P 2 is gm/cc, and P 2 is in atmospheres. According 
to Reference 4-7, the shock layer is frozen if Da^ < 3 x 10 . If 
Da^ » 1, the shock layer is in equilibrium. 

A similar analysis was applied to determine if vibration would 
be in equilibrium behind the normal shock. In this case, Blackman's 
relaxation rate equation (Ref. 48) was used for the excitation of the 
vibrational mode. Thus, 




D V = 


V M 7 2 


-11 1/2 
,1 X 10 t 2 


2 ]r(#) 


where Pg is in atmospheres, Tg is K. 

Inger's results (Ref. 49) were used to calculate the state of 
the boundary layer. The Damkohler number for recombination in the 


boundary layer is 


Cl = Cl a e 


where 


r 2 
2L W T 

Iz ’ p l + 0.47 s?13 ^>*2 + (0.47 S 0 l/3)2 Iz > p 3 


0.4-7 Sc 1 / 5 


Ci . 


Lw = 


30 fcA 

otas/ 


Units and the other parameters, with the exception of I z F _^, I z F ^, 

and l-ZfFj’ have been P reviousl y defined. 811(1 I Z,F 3 

may be found graphically in Reference 49 as a function of T w , T 0 , 
t d> and tD. If Ci < 10“^, the boundary layer may be considered 
frozen (Ref. 45). 

Low-density, second-order viscous effects may also affect the 
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data. Cheng (Ref. 50) presents a figure showing the ratio q/q- 


versus a Reynolds number defined as Ref = 


p t,2 7 H t,l R Jf 




. As the 

^t,2 

Reynolds number decreases, this ratio increases above 1.0 and then 
rapidly decreases as free-molecule flow is approached. This figure may 
be used to obtain a correction to the stagnation point heat -transfer 
rate in the low Reynolds number regime. 

Gas Flow Models 

The problem of the nonequilibrium expansion of gas in a nozzle 
has been investigated for several years. Some of the techniques used 
will be mentioned briefly, but this by no means lists all of the 
contributions. 

Glcrwacki (Ref. 9 ) treated the chemical nonequilibrium behavior 
of oxygen in a dilute mixture of inert N 2 with vibration assumed to be 
in equilibrium. Enkerhus (Ref. 10) treated the approximate nonequilib- 
rium case by allowing vibration and chemistry to freeze at a point 
corresponding to the sudden-freeze point of the oxygen recombination. 
This amounted to an uncoupled sudden-freeze model for chemistry and 
vibration. Again, the composition was O 2 and 0 in a dilute mixture 
oi inert N 2 « 

More exact results have been obtained for more complex gas 
mixtures by the use of electronic computers. Emanuel and Vincenti 
(Ref. 11) used a model for air consisting of N 2 , N, 0 2 , 0, and NO 
which included eight chemical reactions plus the vibrational processes 
of the three diatomic species. Here, again, however, the chemical and 
vibrational modes were uncoupled. 
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Lordi, Mates, and Moselle (Ret. 12) presented a computer program 
for an air model consisting of N 2 , N, 0 2 , 0, NO, N0 + , Ar, and E" with 
11 chemical reactions, although the program can handle 20 species and 
64 reactions. Vibration may be assumed either to be in equilibrium 
with translation and rotation or to be frozen at conditions prior to 
the expansion. 

Tirumalesa (Ref. 13) considered the case of pure oxygen with 
coupled vibration and dissociation. The coupling models were those of 
Treanor (Ref. 51). 

Dressier, French, and Webb (Ref. 14) developed a computer pro- 
gram for the nonequilibrium flow of air, including nonequilibrium 
vibrational and electronic energy relaxation and coupling effects 
between these energy inodes and the chemistry. Rotation, vibration, 
and electronic excitation were considered uncoupled. The program is 
capable of handling 20 species and 50 reactions. 

Thus, it is evident that the available techniques for calcu- 
lating nonequilibrium flows cover a wide spectrum of exactness and 
complexity. An excellent survey of nonequilibrium flow is given in 
Reference 52. This paper notes the status of the problem, both experi- 
mental and theoretical. 

In the present case, advantage was taken of the ready avail- 
ability of the program of Lordi, Mates, and Moselle (Ref. 12) for the 
computation of nonequilibrium chemistry. This program was used in 
conjunction with a sudden-freeze model of N 2 vibration to predict test 
section properties. The procedure will be discussed in detail below. 




Mention will be made of the limiting cases of nonequilfbrium 
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flow in order to present some known bounds to the predicted values of 
flow properties. Thus, the cases of equilibrium and frozen flow will 
be considered briefly since the degree of uncertainty in the chemical 
reaction rates affects the nonequilibrium predictions. 

Equilibrium flow .- The basic assumption in equilibrium flow is 
that all reaction times are zero; that is, the fluid is in complete 
thermodynamic equilibrium at each stage of the expansion. The assump- 
tion of complete thermodynamic equilibrium along with the assumption 
of a constant entropy expansion makes the calculation of equilibrium 
flow simp3e. The specification of one other thermodynamic property 
and the use of a Mollier diagram enables one to solve the equilibrium 
expansion for specified plenum chamber conditions. 

The equilibrium predictions used in the present investigation 
were taken from an existing equilibrium flow expansion computer program 
at Langley Research Center. The species in the simplified equilibrium 
program are N 2 , 0 2 , and 0. Results from this program compared favor- 
ably with the more complete equilibrium solution (more species included) 
of the Cornell Aeronautical Laboratory program (Ref. 12). The simpler 
program was used because entrance to the program was easier since 
normal shock properties were tabulated* 

Nonequilibrium flow .- The nonequilibrium flow solution is much 
more complicated than the equilibrium solution since finite reaction 
rates are included. The assumption of constant entropy of the expan- 
sion can no longer be made and the history of the flow can no longer 
be neglected. 
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The calculation technique of the present study involves using an 

available chemical nonequilibrium program (rotation, translation, and 

vibration in equilibrium) for air (Ref. 12 ) in conjunction with a 

sudden-freeze criterion for the vibrational mode. Ti odel for air 

consists of eight species: N2> N, O2, 0 , NO, Ar, N 0 + , and E". The 

11 reactions considered and the forward reaction rate for each are 

presented in Table I. The program requires a nozzle geometry as 

input. The actual geometry was a circular arc throat (Radius of 

curvature = I.75 inch) followed by a 5° half-angle conical nozzle. 

The actual nozzle was simulated by fitting a polynomial to the equation 

of the circular arc throat in the converging throat region. Downstream 

o 

of the throat, a hyperbolic contour asymptotic to a 5 half-angle cone 
was used. 

The logic involved in computing the nonequilibrium flow is that, 
in most cases, chemistry will freeze before the vibrational mode 
freezes during an expansion. (This is not true in the case of atomic 
nitrogen due to the NO "shuffle" reactions; however, the mass fraction 
of N was two orders of magnitude lower than other species over the 
range of conditions of interest hero. ) Thus, the freezing of both 
chemistry and vibration at the freezing point of vibration provides a 
reasonable estimate of fluid properties (including composition) in the 
te^t section. 




various species in an expansion down the nozzle. The conditions 
chosen represent one of the highest enthalpies encountered during 
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TABLE I. - CHEMICAL REACTIONS USED IN THE NONEQUILIBRIUM 

NOZZLE PLOW PROGRAM 


Reaction 


N 2 + M 2N + M 
N 2 +M£2N+M N 

N 2 + M ^ 2N + M 0 2 ,Ar,0,N0 
0 2 + M 20 + M 0 2 

0 2 + M 20 + M 0 

0 2 + M £ 20 + M N 2 .Ar,N,N0 


®act» 


mole-K 


-1.5 -2.21*99 x io5 

1.5 x 10 cc - 1.5 -2.2499 x 105 
9.9 x 10 20 - 1.5 -2.2499 x io5 

3.6 x 10 21 - 1.5 - 1.1796 x 105 

2.1 x 10 1 ® - 0.5 - 1.1796 x 10 5 

1.2 x 10 21 - 1.5 - 1.1796 x 105 


N0+M??N + 0+ M N 2 ,0 2 ,Ar,N,0,N0 5.2 x 10 21 - 1.5 -I .4996 x 105 
0 2 + N ^ 0 + NO — 1.0 x 10 12 0.5 -6.2 x 103 

N g + 0 2N0 + N — 5 x lO 1 ^ 0 -7.552 x 10^ 

Ng + 0 2 2 NO + NO — 9.1 x lO 2 ^ - 2.5 -1.2912 x 105 


N 0 + + E" £ N + 0 


1.8 x 10 21 l-1.5 


3 AiT - exp|^± 


T = °K 


\ = 1.98647 - 

5 , m 

= car /mole-sec 


mole-°K 
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the tests. Ifte solid lines are for chemical nonequilibrium with other 
modes in equilibrium, while the dashed lines represent complete thermo- 
dynamic equilibrium. Both sets of data were taken from the computer 
program described i** Reference 12. The vertical dashed line indicates 
the area ratio where vibration is predicted to freeze according to the 
criterion to be described below. Note that all species except N, N0 + , 
and E" are frozen at this point. The nonfrozen species are present 
only in very small quantities. 

The criterion used to establish the vibrational freezing point 
was that of Fhinney (Refs. 55 and 54). The freezing point is defined 
as the location where 


dE 

dt 


E 


( 25 ) 


E is the equilibrium vibrational energy of the gas, and t is the 
vibrational relaxation rate. This criterion is similar to that for 
recombination freezing given in Reference 55* dE/dt was approximated 
as 


AE 

E 3 “ Ei 


(Ej 

At 

tj - tx 





x 5 - X 1 


( 26 ) 


where E and t v were calculated at point 2, 1 was upstream, and 

5 downstream of 2 (the point in question). 

The relaxation rate used was that for N 2 in air. This rate was 
determined by using measurements of stained during the present 

studies and fitting various rate equations of the form 


1 
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* = Al exp (^j) (27) 

where T = °K, P = atm, to the experimental data by use of the N 2 
vibrational nonequilibrium program of Erickson (Ref. 56 ). The vibra- 
tional temperature data will be discussed under results and discussion. 

The vibrational modes of the diatomic molecules N 2 , 0 2 , NO, and 
N0 + do not necessarily relax at the same rate. However, due to uncer- 
tainties in the relaxation rates of the other diatomic species, as well 
as the fact that N 2 forms the bulk of the diatomic molecules present, 
the use of the N 2 vibrational rate for the gas seems justified. 

Once the freezing criterion set forth above was satisfied, the 
chemical nonequilibrium program was terminated and all fluid properties 
were tabulated at that station. From this point on, the chemical com- 
position and the vibrational energy were considered frozen. 

This assumption of frozen flow allows the use of perfect gas 
relations with an effective ratio of specific heats to calculate the 
flow properties downstream of the freezing point. This effective ratio 
of specific heats is dependent on the composition at the freezing point 
and may be defined as 


where 


and 


5 ^ + r a 

7eff ^ 

(28) 

* R vK 2 + + ’'NO + t 'noO 

(29) 

R A “ R uQ'e“ + ’'n + r| 6 + n Ar3 

(JO) 
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The reference properties used in the frozen expansion are those at the 
freezing station. Properties at this station axe denoted by a super- 
script (')• The perfect gas property relationships may be found in 
•several references (e.g., Refs. 8 and 57 ). 

The properties at different points downstream of the freezing 
point are calculated by choosing a temperature T and calculating the 
other properties which are needed for comparison with the experimental 
measurements. Thus, in terms of T and the superscript (’) freezing 
point properties, expressions for P, p, V, A/A*, and M are as 
follows : 


P / T VeffAeff 1 


1 . (15 

P' \T 7 


P* IT 1 


TWreff" 1 


v = J(v* T + 2 c P eff ( T * - T) 


■ (ifeMW 


7 eff( R A + R m) T 


For all of the present tests, the shock layer around the pitot 


probe was chemically frozen since the DamkShler number for dissociation 
(Eq. (20)) was less than 3 x 10“ 2 . Vibration was also frozen according 
to the Damkohler number for vibration (Eq. (21)). Thus, chemistry and 
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vibration were assumed constant at free-stream values in the shock 
layer and the same effective ratio of specific heats was used behind 
the shock as in the free stream. So, the flow across the normal shock 
was a perfect gas flow with an effective ratio of specific heats, 

7 e ff , Equation (28). Expressions for the properties behind the shock 
and at the stagnation point may be found in Reference 57* 

Frozen flow .- In frozen flow, as considered here, chemical reac- 
tion and vibrational relaxation times are infinite. Consequently, for 
chemically and vibrationally frozen flow, the composition of the flow 
as well as the chemical and vibrational energy remain constant and do 
not contribute to the expansion. Again, the flow expansion may be 
calculated as that of a perfect gas with an effective ratio of specific 
heats. 

In the present studies, the limiting case of frozen flow was 
taken to be that with chemistry and vibration frozen at the values in 
the plenum chamber. This case along with the equilibrium limiting case 
indicates the bounds of any possible flow property variation. 

Once again, normal shock properties were calculated by perfect 
gas formulae using an effective ratio of specific heats since the flow 
was frozen behind the shock. No equations for the normal shock 
properties or the free-stream properties will be given since they are 
all found in Reference 57* 




CHAPTER V 

RESULTS AND DISCUSSION 

The experimental and theoretical data of these studies are 
presented in both dimensional and dimensionless form with the excep- 
tion of the calorimeter data, the stagnation point heat-transfer data, 
and the axial pitot pressure distributions. Stagnation chamber 
properties are used to nondimensionalize the test stream properties. 
Total enthalpy is nondimensionalized by R S T 0 , while free-stream 
velocity is nondimensionalized by the limiting velocity for complete 
conversion of thermal energy to kinetic energy. The "(a)" part of 
a figure will always be the dimensional data, while the "(b)" part 
will be the dimensionless data. 

Calorimeter Data 

A knowledge of the total enthalpy is essential in a study of 
this type where one is attempting to determine the degree of nonequilib- 
rium of the flow. A theoretical technique such as the sonic-throat 
method must be experimentally verified, since this method is dependent 
on the state of the gas as well as the throat discharge coefficient. 

In addition, the geometry of the heater affects the prediction of this 
method in swirling flows. Proof of this has been obtained by the 
author by simulating such swirling flows with room-temperature air 
(to be published). 
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The enthalpy measured by the total calorimeter is compared to 
the equilibrium sonic-throat method in Figure 15. These data were 
taken with two different outer electrodes as noted. Total calorimeter 
enthalpy ranges from 1,180 Btu/lbm to 4,020 Btu/lbm. 

The comparison of H t#sonic throat and H^ CAL closely 
parallels the line of perfect agreement to values of 2,500 Btu/lbm. 

At higher enthalpies, the sonic-throat method increasingly overpredicts 
the enthalpy. This may be caused by increasing swirl in the flow due to 
rotation of the arc as enthalpy increases (arc speed remained constant 
for each electrode). The evidence of swirl or vortex flow was collected 
by measuring the pressure on the center of the flow (on the inner 
electrode) and on the wall of the contamination trap (4 inches down- 
stream from the tip of the inner electrode). As enthalpy increased, 
the wall pressure grew increasingly larger than the centerline pressure 
(see Fig. 16). Since wall pressure was used in the sonic-throat tech- 
nique, such a trend would give the results shown in Figure 15. 

The reason for the increasing swirl is not definitely known, 
but it could be due to the effects of increasing flow contamination 
with increasing enthalpy. Such a trend has been measured in other 
facilities (Ref. 58). Attempts were made to measure the amount of 
contamination as a function of enthalpy during the present studies. 

Scatter of these data was severe, although the trend was one of 
increasing amounts of contamination as total enthalpy increased. 

Disagreement of the sonic -throat prediction with the calorimeter 
at the lower enthalpy levels where swirl did not seem to be a factor 
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Variation of the difference between arc chamber 
outside wall pressure with total enthalpy* 
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could be due either to nonequilibrium effects or to effects of arc 
radiation on heat losses to the throat. Experimental data obtained 
from the calorimeter tests are tabulated in Table II of the Appendix. 

Stagnation Point Heat Transfer 

Stagnation point heat-transfer rate is quite often used to pre- 
dict total enthalpy in the test section of wind tunnels. However, as 
mentioned previously, this method is rather uncertain in low-density 
flows where chemistry and degree of catalycity of the gauge surface 
must be considered. This situation existed in the present studies 
as well as some uncertainty of the effective radius of the hemispherical 
nose section. 

The stagnation point heat -transfer data are shown in Figure 17. 
The data are plotted as q/ v/Pt,2' n0 nose radius has been 

introduced into the experimental values. Zoby's theory for a com- 
pletely catalytic wall with no consideration of second-order effects 
is shown for the limits of the effective nose radius. The calculations 
of the Damkohler number for recombination in the boundary layer indi- 
cated that the boundary layer was frozen (c^ < 10“ i ). Thus, the 
uncertainty of the speed of reaction at the gauge surface must be 
considered in comparing the heat-transfer data with theory. This 
uncertainty was heightened in the present tests by deposits of contam- 
inants on the surface. At the highest enthalpy of these tests, a 
spread of 30 per cent would theoretically exist between a completely 




catalytic and a completely noncatalytie surface (by Eq. (14)). The 
present data should fall within these limits, and they do. 
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Previous analyses of contamination deposits in tests of this 
type indicated oxides of copper plated on the surface. Therefore, 
values of were calculated for cupric oxide, using the information 
on surface recombination coefficients tabulated in Reference 43. At 
the highest value of enthalpy of the present tests, the value of heat- 
transfer rate for a surface covered with cupric oxide lies almost 
exactly between the values for a perfectly catalytic and a perfectly 
noncatalytie surface, or 15 per cent higher than the perfectly non- 
catalytic surface. This percentage decreases with decreasing enthalpy 
until essentially no difference is indicated between a catalytic and 
a noncatalytie surface at an enthalpy of approximately 1,500 Btu/lbm. 

According to the theory of Cheng (Ref. $ 0 ), all of the present 
data are approximately 12 per cent above the high Reynolds number 
boundary-layer prediction because of the second-order viscous effect. 
The scatter of the experimental data that have been compared with this 
theory makes the actual magnitude of the correction somewhat 
questionable. 

A more direct comparison with the catalytic, high Reynolds 
number heating-rate theory could be made if the data were corrected 
for chemical and second-order viscous effects. This would involve 
correcting the experimental data shown in Figure 17 by approximately 
-12 per cent to account for the viscous effect. No correction for 
chemistry and wall catalytic effects would be necessary below 
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approximately 1,500 Btu/lbm. Above this level, the data should be 
correct ew by an increasing percentage for surface catalytic effects 
(+15 per cent at the highest enthalpy). This is assuming that the 
surface is plated with cupric oxide. Thus, the data would show a net 
change ranging from a 12 per cent decrease at the lowest enthalpy to a 
5 per cent increase at the highest enthalpy. These proposed corrections 
are only noted and are not applied to the data. 

As can be seen from Figure 17, the limits of uncertainty on the 
nose radius of the heat-transfer probe impose about a 20 per cent 
spread on total enthalpy at a given value of heat-transfer rate. Since 
the data (either corrected or uncorrected for chemical and second-order 
viscous effects) fall between these limits, the value of total enthalpy 
in the free stream is certainly confirmed much closer than this 
(probably within ±15 per cent). 

These tests emphasize the problems associated with accurately 
measuring heat-transfer data in low-density, high-enthalpy facilities. 
Also, the critical importance of knowing the exact nose shape is 
obvious. Experimental data obtained from the stagnation point heat- 
transfer tests are tabulated in Table III of the Appendix. 

Vibrational Temperature 

The vibrational temperatures were determined by comparing the 
relative intensities of the (0,1) and (1,2) vibrational bands of Ng. 

The vibrational temperatures of N2 measured in air are shown in 
Figure l8(a) as a function of T^.^ and in Figure 18(b) as a function 


m 
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Figure l8.- Measurements of N 2 vibrational temperature illustrating the 

nonequilibrium nozzle flow. 
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of Ej. -,/R s T 0 * The fact that the vibrational temperature is not in 
equilibrium is shown, and these data are used to establish a rate 
equation for the vibrational relaxation of N2 in air as discussed in 
the section on nonequilibrium flow theory. This equation was used to 
determine the effect of vibrational freezing on the other gas dynamic 
parameters of the flow. The rate equation that best fit the experi- 


mental data was 


T y p * l.l x 10' 12 exp 


fe$j| 


atm, sec 


( 36 ) 


with T in °K. These results on the N2 vibrational relaxation rate 
are in agreement with the experimental results of References 7 and 18 
that the relaxation rate is faster in an expanding flow than behind a 
normal shock (Ref. 48 ). The vibrational temperature data obtained 
from these tests are tabulated in Table IV of the Appendix. 


Rotational Temperature 

The experimental values of rotational temperature are compared 
with theory in Figures 19(a) and 19 (b). Although the ordinate is 
labeled the rotational temperature of it is understood that rota- 
tion and translation v are represented by the same kinetic temperature 
and that this temperature is the static temperature of the gas mixture. 
The curves shown are the predictions of the theories of complete 
thermodynamic equilibrium, chemical nonequilibrium with vibration in 
equilibrium, uncoupled chemical and approximate vibrational nonequilib- 
rium, and frozen chemistry and vibration. 


(a) Dimensional data. 


Figure 19 •- A comparison of measured Ng rotational temperature with 
theoretical predictions of the nozzle expansions of gas models 
ranging from frozen to equilibrium. 
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Two sets of symbols are shown on the figure. The open symbols 
are the rotational data temperature obtained using the Muntz model 
(Ref. 2k) and the solid symbols are the data reduced according to the 
model proposed by Hickman (Ref. 29). Temperatures obtained using the 
Hickman model are in better agreement with the chemical and vibrational 
nonequilibrium theory, although the Muntz model is more accepted and 
has received wider usage. As noted earlier, the Muntz model over- 
estimates the temperature, especially at low temperatures, although it 
becomes increasingly accurate as temperature increases. The per cent 
difference between the two models is 20 per cent at 100° R and 5 per 
cent at 400° R. 

The Hickman model is considered to be the more accurate of the 
two models in lieu of its performance at the lower temperatures and in 
view of the results of Ashkenas' work in this density range. Note 
that at the lower temperatures, the data as resolved by the Muntz model 
lie either above or nearly coincident with the equilibrium theory up 
to a total enthalpy of approximately 1,500 Btu/lbm. Consideration of 
only the Muntz model could lead to the conclusion that vibrational 
equilibrium is unimportant. However, from the present calculations 
and the work of Erickson (Ref. 5 6 ), the degree of vibrational non- 
equilibrium must clearly be taken into account. Hence, in speaking of 
the rotational temperature data henceforth, the reference will be to 
the data as reduced by Hickman's model. 

Generally, the data lie slightly above the chemical and vibra- 
tional nonequilibrium theory over the entire range of total enthalpy. 
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It is difficult to comment on the model used to simulate vibrational 
nonequilibrium, since the data lie equally close to the curve for 
chemical nonequilibrium with vibration in equilibrium. 

One conclusion that may be drawn from this figure is that the 
gas is in a nonequilibrium state. The rotational temperature reduction 
methods are certainly sufficiently accurate to confirm that the flow 
is neither in equilibrium nor completely frozen. Also, within the 
probable accuracy of the data reduction models, the chemical and 
vibrational nonequilibrium theory proposed in this report is' in reason- 
able agreement with the data. Final judgment of this theory can only 
be made after comparison with the other measured parameters. The 
rotational temperature data are tabulated in Table V of the Appendix. 

Nozzle Wall Static Pressures 

The nozzle wall static pressure data are presented in 
Figures 20(a) and 20(b) as a function of total enthalpy. The open 
symbols were obtained with the l/4-inch-diameter flat-faced cylindrical 
pitot probe located on the tunnel centerline, while the flagged-open 
symbols were made with no probe in the tunnel. Thus, it is evident 
that the presence of the probe lid not affect the wall static 
pressures. 

The data have been corrected for the heat-transfer orifice 
effect (Refs. 33 and 34). A3 per cent to 12 per cent upward correc- 
tion with increasing enthalpy was necessary. No corrections were made 
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(a) Dimensional data* 


Figure 20.- A comparison of measured wall static pressure with 
theoretical predictions of the nozzle expansions of gas 
models ranging from frozen to equilibrium. 






Figure 20.- Concluded 





for static pressure variation across the flow core as was done in 
Reference 7, since no clear-cut method was available to make such 
a correction. 

The need for some further correction to the pressure data is 
obvious if one considers the lower enthalpy data. Note that the spread 
between the frozen and equilibrium theories is approximately 13 per cent 
at 1,000 Btu/lbm, while the experimental data lie about 3° per cent 
above the equilibrium theory. The experimental measurements should be 
accurate to within approximately 5 per cent at this pressure level. 

Thus, neither the accuracy of the data nor the possible spread due to 
the state of the gas can account for the observed disagreement with 
theory. 

The data lie considerably above nonequilibrium theory over the 
entire range of enthalpy and cross equilibrium theory at approximately 
1,900 Btu/lbm. Compression waves caused by the contouring effect of 
the boundary layer and/or lateral flow gradients could tend to make the 
wall pressure higher than centerline pressure. 

Thus, the wall static pressures are inadequate (at this high 
area ratio) for the determination of the proper expansion model. A 
"calculated" experimental centerline static pressure using the measured 
rotational tenperature and the measured free- stream density from the 
mass flow probe (to be presented in the next section) could be calcu- 
lated from the thermal equation of state, P * ZpR fl T, since only 

varies from 1.00 to 1.14. However, the uncertainty involved with the 
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rotational temperatures negates this approach. The experimental static 
pressure data are tabulated in Table VI of the Appendix. 

Static Density 

The free-stream static density, as determined from the mass 
flow probe, is shown as a function of total enthalpy in Figures 21(a) 
and 21(b). The chemical and vibrational nonequilibrium theory agrees 
well with the experimental values of density above approximately 
2,000 Btu/lbm. Below 2,000 Btu/lbm, the data gradually trend upward 
to the frozen expansion theory. The spread between nonequilibrium 
and frozen theory is on the order of 10 per cent in this enthalpy 
range, however, and this upward trend could be simply data scatter. 
Another possibility for the discrepancy is the factor B used in the 
relationship between pitot pressure and the product of free-stream 
velocity and pressure, that is, ^t,2 = In these tests, B 

was assumed to be constant at O.925, although it is weakly dependent 
on the effective value of the ratio of specific heats. 

The density data generally confirm the nonequilibrium theory 
proposed in this report. Although the lower enthalpy data are diffi- 
cult to interpret, the spread between frozen and equilibrium flow at 
the higher enthalpies and the good agreement of the experimental data 
with the chemical and vibrational nonequilibrium theory in this range 
enable one to arrive at the preceding conclusion. The experimental 
density data are tabulated in Table VII of the Appendix. 
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Velocity 

The velocity measurements obtained from the mass flow probe are 
shown in Figures 22(a) and 22(b) as a function of total enthalpy. The 
chemical and vibrational, nonequilibrium theory generally is in much 
closer agreement with the experimental data than are the other theoret- 
ical gas models. At the lower enthalpies, as in the case of density, 
the velocity data fall onto the frozen theory prediction. The possible 
source of this behavior has been discussed previously in the section 
concerning free-stream density. 

The good agreement with nonequilibrium theory serves to give 
one confidence in the present theoretical formulation. The allowance 
for vibrational freezing is reasonably verified when one compares the 
chemical nonequilibrium theory (vibration in equilibrium) with the 
experimental data. The experimental velocity data are tabulated in 
Table VII of the Appendix. 

Pitot Pressure Gradients 

All of the pitot pressure data accumulated during the present 
experiments are shown in Figure 23 as a function of total enthalpy. 
These data were obtained with the various probes noted on the figure 
at locations of -6, 2-1/1+, and 13 inches from the nozzle-test section 
juncture. 

The composite of these data is effectively an axial survey of 
the pitot pressure on the centerline of the tunnel. With a knowledge 
of the proper gas expansion model, one can specify flow gradients, such 
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as the longitudinal Mach number gradient. The Mach number gradient 
is approximately 0.06 per inch at 7 00 Btu/lbm and approximately 
O.O 35 per inch at 3 >500 Btu/lbm. 


CHAPTER VI 


CONCLUDING REMARKS 

A series of tests have been performed to calibrate and diagnose 
the flow in an electric arc heated test stream of air. Tests were made 
over a range of plenum chamber conditions with plenum chamber enthalpy- 
varying from 700 to 3*510 Btu/Ibm and plenum chamber pressure varying 
from 4.20 to 9*38 atmospheres. Parameters measured during the tests 
were N^ vibrational temperature, rotational temperature, wall 
static pressure, static density, velocity, total enthalpy, and hemisphere 
stagnation point heat-transfer rate. 

Based on the results of these tests, the following conclusions 
can be made: 

(1) The rotational temperature data indicate that the flow is 
in a nonequilibrium state since the data fall between the prediction 
of the chemical nonequilibrium theory with and without vibrational 
nonequilibrium. Results frcm Hickman's data reduction technique for 
obtaining the rotational temperature were in closer agreement with the 
chemical nonequilibrium sudden-vibrational freezing theory than were 
the results of Muntz's method. The percent difference between the two 
methods ranged frcm 20 percent at 100° R to 5 percent at 400° R. 

(2) Measured density and velocity (obtained frcm the mass flow 
probe) generally agreed with the predictions of the chemical nonequi- 
librium, sudden- vibrational freezing theory. The consistency of the 
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experimental data and the good agreement of these data with theory 
support the design of the mass flow probe with equal interior and 
exterior tip angles. 

(5) Measurements of N 2 vibrational temperature in air indicated 
vibrational freezing just downstream of the nozzle throat. The data 
were used to obtain a vibrational relaxation rate equation. The 
vibrational temperature data also supported the results of other inves- 
tigators that the vibrational mode relaxes faster during an expansion 
than behind a normal shock. 

(k) Wall static pressures were inconclusive in establishing 
the proper gas flow model to use for the nozzle expansion. Weak 
compression waves caused by boundary layer countouring may be the cause 
of the high measured pressure data. However, the static pressure data 
definitely confirm that the flow is in a nonequilibrium state. 

(5) The sonic-throat technique for predicting stagnation 
enthalpy should be experimentally verified before relying upon the 
method. This is evident from the present tests where such factors as 
swirl, heater geometry, and the thermodynamic state of the gas cause 
deviation from equilibrium sonic-throat predictions. 

(6) The comparison of stagnation point heat-transfer data with 
theory was affected by uncertainty of the surface catalytic efficiency, 
the magnitude of the second-order viscous effect, and the uncertainty 
of the hemisphere nose radius. However, the free- stream total enthalpy 
could still be verified to approximately tl5 per cent by the use of the 


heat -transfer measurements. 
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Thus, the conclusions generally point to acceptable agreement of 
the various measured experimented data with the predictions of a chem- 
ical nonequilibrium, vibrational sudden-freeze technique for calculating 
the nozzle expansion of the air. Several pertinent free- stream param- 
eters can be determined by using this flow expansion gas model. At a 
location 2-1/4 inches from the nozzle test section juncture, the Mach 
number varies from 11.2 to 12.3 over a total enthalpy range from 700 
to 3,500 Btu/lbm. This corresponds to a stagnation chamber pressure 
range from 4.20 to 9*40 atmospheres. Reynolds number per foot at this 
location varies from 3 * 10^ to 1 x lo\ while the Mach number axial 
gradient varies from 0.06 per inch to 0 035 per inch over this same 
range of supply conditions. 
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APPENDIX 
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TABLE III.- EXPERIMENTAL RESULTS FROM THE STAGNATION POINT 

HEAT-TRANSFER RATE TESTS 


Test no. 

p t,i 

(atm) 

«t,i 

(Btu/lbm) 

T t,i 

<°R) 

• 

q 

(Btu/ft 2 -sec) 

P t,2 

(atm) 

15 

4.27 

840 

3130 

11*9 

6.28 x 10"3 

16 

4.67 

990 

3635 

14.6 

6.67 x io ~3 

17 

5.44 

1225 

4365 

18.1 

6.86 x 10“3 

18 

5.27 

1250 

4430 

21.9 

7.78 x io -3 

19 

5.78 

1370 

4750 

23.I 

8.01 x 10-3 

20 

VJ1 

CD 

VJ1 

1500 

5075 

28.9 

8.89 x 10-3 

21 

6.04 

1510 

5H0 

27.O 

8.62 x 10"3 

22 

6.31 

1600 

5330 

28.7 

8.88 x 10"3 

23 

6.65 

1680 

5510 

27.I 

8.91 x 10“3 

24 

6.6 5 

1760 

5690 

29.9 

9.31 x 10"3 

25 

6.56 

1790 

5740 

35*0 

9.78 x 10-3 

26 

6.68 

1890 

5940 

36.9 

9*78 x 10'3 

27 

903 

2420 

6770 

51.6 

H.3 x 10”3 

28 

9.98 

2680 

7130 

58.4 

11.5 x io "3 

29 

10.08 

2805 

7280 

62.4 

11.9 x io -3 

30 

7-33 

2855 

7255 

58.0 

9.89 x io "3 

31 

10.08 

2900 

7390 

71*3 

13.2 X 10 -5 

32 

9.31 

2940 

7435 

73-6 

12.2 x io "3 

33 

4.60 

2970 

7290 

52.0 

7.21 x io "3 










Test no 

34 

35 

36 

37 

38 

39 

40 

41 

42 

^3 

44 


im 

H t,i 

(Btu/lbm) 

4.35 

700 

4.76 

905 

5.51 

1200 

5.66 

1280 

6. 02 

1410 

6.96 

1650 

6.61 

1855 

8.46 

1970 

7*69 

2240 

8.15 

2250 

8.25 

2520 

8.46 

2520 

9*03 

2980 


45 

46 

47 


9.51 


3130 








TABLE V. - EXPERIMENTAL RESULTS FROM THE N 2 ROTATIONAL 

TEMPERATURE TESTS 


Test no. 

p t,i 

(atm) 

H t,l 

(Btu/lbm) 

T t,l 

(°R) 

Tr,Hp,« 

(M 

P t,2 

(atm) 

48 

4.20 

700 

2665 

92.9 

5.10 x io “5 

49 

4.86 

940 

3475 

133.2 

5.56 x 10“3 

50 

5*17 

1080 

3940 

116.7 

6.21 x 10’5 

51 

6.15 

1360 

4735 

173-0 

6.90 x 10'5 

52 

6.89 

1380 

4770 

180.9 

7.19 x 10“2 

53 

6.63 

1440 

4930 

192.2 

7.42 x 10-5 

54 

8.42 

2380 

6695 

284.1 

9.10 x 10‘5 

55 

8.49 

2450 

6803 

273.9 

9.63 x io ”5 

56 

8.76 

298O 

7470 

312.2 

9.87 x 10"5 

57 

9.78 

3220 

7810 

385.6 

IO.91 x 10"5 

58 

9.72 

3310 

7920 

305.8 

IO.89 x 10 '5 









Test no. 

ess 

59 

4.13 

*60 

4.65 

61 

5*56 

*62 

6.54 

63 

6.75 

64 

6.99 

65 

8.29 

*66 

8.92 

67 

8.98 

68 

9.10 

69 

9-52 


H t,i 

(Btu/lbm) 

T t,i 

(°R) 

780 

2990 

1010 

5710 

1410 

4860 

1830 

5815 

1900 

5940 

1920 

5975 

2870 

73^5 

3070 

7595 

3080 

7615 

3080 

7615 

3115 

7670 


(atm) 

6.11 x 10*5 
7.26 x 10*5 

7.95 x 10*5 
9.22 x 10*5 
8.80 x 10*5 
9.25 x 10“5 

10.6 x 10*5 

11.6 x 10’5 

11.8 x 10*5 

11.3 x 10*5 
11.5 x 10*5 



7.15 x 10*5 


9.43 x 10*5 

11.3 x 10*5 

11.4 x 10*5 
13*3 x 10*5 

14.3 x 10*5 

13.8 x 10*5 

13.8 x 10*5 
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TABLE VII.- EXPERIMENTAL RESULTS FROM THE MASS FLOW 

PROBE TESTS 


Test no. 

mam 

1 

H t,l 

(Btu/lbm 

■pg 

■n 

Poo 

( slugs/f t3 ) 

Vco 

(ft/sec) 

imam 
1 1 

TO 

5-54 

1290 

4555 

2.83 x 10*7 

7280 

6.54 X 10*5 

71 

6.22 

1610 

5345 

2.78 x io*7 

7810 

7.41 x io*3 

72 

6.51 

1725 

5580 

2.79 x io"7 

786O 

7-54 x 10*3 

73 

7.27 

1730 

5580 

2.73 x 10*7 

8290 

8.21 x 10*3 

74 

7-35 

2080 

6175 

2.45 x 10*7 

9040 

8.73 x 10*3 

75 

7.^0 

2080 

6175 

2.47 x io*7 

8990 

8.72 x 10*3 

76 

8.04 

2460 

6805 

2.34 x 10*7 

9570 

9*36 x io*3 

77 

8.45 

2520 

6875 

2.46 x 10*7 

9550 

9.80 x io*3 

78 

8.08 

2560 

6 930 

2.36 x io*7 

9760 

9.81 x 10*3 

79 

8.26 

2720 

7165 

2.26 x io*7 

9960 

9-76 x io*3 

80 

10.05 

3290 

7900 

2.06 x io*7 

11430 

11.8 x 10" 5 

81 

9.28 

3510 

7940 

2.08 x io"7 

11300 

11.6 x 10-3 

82 

9.58 

3510 

8225 

2.15 x io“7 

11280 

12.0 x io*3 






